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The bulk viscosity in spin-one color-superconducting strange quark matter is calculated by taking 
into account the interplay between the nonleptonic and semi-leptonic week processes. In agreement 
with previous studies, it is found that the inclusion of the semi-leptonic processes may result in 
non- negligible corrections to the bulk viscosity in a narrow window of temperatures. The effect is 
generally more pronounced for pulsars with longer periods. Compared to the normal phase, however, 
this effect due to the semi-leptonic processes is less pronounced in spin-one color superconductors. 
Assuming that the critical temperature of the phase transition is much larger than 40 keV, the main 
effect of spin-one color superconductivity in a wide range of temperatures is an overall increase of 
the bulk viscosity with respect to the normal phase. The corresponding enhancement factor reaches 
up to about 9 in the polar and yl-phases, about 25 in the planar phase and about 29 in the CSL 
phase. This factor is determined by the suppression of the nonleptonic rate in color-superconducting 
matter and, therefore, may be even larger if all quark quasiparticles happen to be gapped. 
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I. INTRODUCTION 

Understanding the physical properties of baryonic 
matter above nuclear saturation density is one of the 
fundamental challenges in modern nuclear astrophysics. 
Many aspects of neutron stars (e.g., the mass-radius re- 
lation, cooling and rotational dynamics, glitches and pul- 
sar kicks) depend on these properties. For example, the 
equation of state of supranuclear baryonic matter plays 
the key role in determining the maximum possible mass 
of neutron stars. The harder (softer) equation of state is, 
the larger (smaller) maximum mass can be. The equation 
of state of dense baryonic matter is also one of the essen- 
tial ingredients that determines the dynamics of core- 
collapse supcrnovae and, in turn, the mass distribution 
of black holes in the Universe [l[ . Several characteristics 
of the gravitational wave emission from merging neutron 
stars [2[, which may be observed, e.g., by advanced Laser 
Interferometer Gravitational Wave Observatory (LIGO), 
are sensitive to the details of the equation of state. 

From the theoretical viewpoint, currently there is no 
consensus even regarding the qualitative state of mat- 
ter at the highest densities reached in stellar cores. The 
most conservative possibility is that such matter is made 
of only nucleonic degrees of freedom. The study of several 
neutron stars in Ref. Q, for example, docs not exclude 
such a possibility, although a phase transition may be 
in agreement with their analysis, provided no extreme 
softening of the equation of state occurs. Another astro- 
physical determination of masses and radii of three neu- 
tron stars in Ref. Q suggests, however, that the actual 
equation of state is too soft to be pure nucleonic. Such 
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contradictory interpretations are representative and show 
that the current knowledge is too limited to settle the is- 
sue. Theoretically, it may be also appropriate to mention 
that the observables associated with the equation of state 
alone have limited power to probe the actual nature of 
dense matter [1, Q. In fact, a true insight regarding the 
stellar interior may require a comprehensive understand- 
ing not only of the thermodynamical, but also transport 
properties and neutrino emission rates of various possible 
states of supcrdense matter. 

In this paper, we assume that baryonic matter at the 
highest stellar densities is deconfincd quark matter. The 
possible formation of quark matter in stars is an old hy- 
pothesis that dates back to the time when the concept 
of quarks was first introduced [8, 9]. This is also sup- 
ported by general considerations jlOj] based on the prop- 
erty of asymptotic freedom in quantum chromodynam- 
ics (QCD) [TTj . The main uncertainties of this scenario 
are (i) the value of the critical density, at which the de- 
confinement transition occurs, and (ii) the actual highest 
density reached in stars. If quark matter is formed, as we 
assume here, it is also likely to be a color-superconductor 
EH- (For reviews on color superconductivity and 
its general effects on stellar properties see, for example, 
Refs. [T^4l9| .) In this paper, in particular, we concen- 
trate on the scenario, in which color superconductivity is 
due to same-flavor, spin-one Cooper pairing [2fj| - |23j . 

The fact of liberation of quark degrees of freedom and 
the formation of a color-superconducting state of mat- 
ter is likely to be revealed through a detailed study of 
the observational features of neutron stars. For exam- 
ple, one promising class of observables is related to the 
rates of weak processes. Such processes are known to 
be responsible for the cooling rates 1241 and damping of 
the rotational (r-mode) instabilities [25( in stars. The 
cooling is primarily determined by the neutrino emission 
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rate, while the damping of r-modes is controlled by the 
viscosity of dense matter [2(| . 

The bulk viscosity in the normal phase of three-flavor 
quark matter is usually dominated by the nonleptonic 
weak processes [27l - l33j . It was argued in Ref . [13] , how- 
ever, that the interplay between the semi-leptonic and 
nonleptonic processes may be rather involved even in the 
normal phase of quark matter. Indeed, because of the 
resonance-like dynamics responsible for the bulk viscos- 
ity and because of a subtle interference between the two 
types of the weak processes, a larger rate of the nonlep- 
tonic processes may not automatically mean its dominant 
role. In fact, it was shown that the contributions of the 
two types of weak processes are generally not separable 
and that, for a range of parameters, taking into account 
the semi-leptonic processes may substantially modify the 
nonleptonic result [34j |. 



In this paper, we extend the analysis of Ref. 3J| and 
study the effect that spin-one color superconductivity has 
on the bulk viscosity when the interplay between the two 
types of weak processes is carefully taken into account. 
(For calculation of the bulk viscosity in other color super- 
conducting phases see Refs. |35l - [39| .) The necessary in- 
gredients for the calculation of the bulk viscosity are the 
rates of semi-leptonic (Urea) and nonleptonic weak pro- 
cesses. While the needed rates for the semi-leptonic pro- 
cesses in several spin-one color superconducting phases 
were obtained several years ago in Ref. H3, the corre- 
sponding rates of the nonleptonic processes remained un- 
known until very recently (4p| . Here we utilize both to 
obtain the bulk viscosity. 

The rest of the paper is organized as follows. The gen- 
eral formalism for the calculation of the bulk viscosity in 
strange quark matter with several active weak processes 
is reviewed in Sec. [TTJ This formalism is then used in 
Sec. IIIII to obtain our main results for the bulk viscosity 
as a function of temperature and the frequency of den- 
sity oscillations. There we also study the enhancement 
effect of color superconductivity on the bulk viscosity, 
as well as the interplay of semi-leptonic and nonleptonic 
processes. In Sec. IIV1 we discuss the results and their 
potential implications for the physics of compact stars. 
Two appendices at the end of the paper contain our fits 
for the numerical suppression factors of the semi-leptonic 
and nonleptonic rates in spin-one color superconducting 
strange quark matter. 



II. FORMALISM 

In this study, in order to calculate the bulk viscosity 
in the presence of several types of active weak processes, 
we follow the general formalism of Ref. [34| . We assume 
that small oscillations of the quark matter density are 
described by 5n = Sn Re(e tut ) where 5n is the magni- 
tude of the oscillations. For such a periodic process, the 
bulk viscosity £ is defined as the coefficient in the expres- 
sion for the energy-density dissipation averaged over one 



period, t = 2-k /u, 

(£diss) 



dt (V • vf , 



(1) 



where v is the hydrodynamic velocity associated with 
the density oscillations. By making use of the continuity 
equation, n + n V • v = 0, we derive 



,c \ C^ 2 / Sn 

(fdiss)--— I — 



(2) 



Such an energy-density dissipation of a pulsating hydro- 
dynamic flow is the outcome of a net work done on a 
macroscopic volume over a period of the oscillation, 



<4 



PVdt, 



(3) 



where V = 1/n is the specific volume. By matching the 
hydrodynamic definition in Eq. ([2]) with the relation in 
Eq. ([3]), we derive the expression for the bulk viscosity, 



(—)' 
\6n J 



PVdt. 



(4) 



The dominant mechanism behind the bulk viscosity is re- 
lated to weak processes [2614291 ] . A periodic oscillation of 
the density is responsible for an instantaneous departure 
from /3-equilibrium in the system. As a result, the for- 
ward and backward weak processes (e.g., u + d s + u 
and s + u — > u + d), which have equal rates in equilib- 
rium, become unbalanced. Their net effect is to restore 
the equilibrium composition. However, since the weak 
rates are relatively slow, a substantial time lag between 
the oscillations of the fcrmion number density (and, thus, 
the specific volume) and the chemical composition (and, 
thus, the pressure) develops. If the resulting relative 
phase shift of the two oscillations is A(f>, one finds from 
Eqs. ([3]) and (|4]) that the corresponding energy dissipa- 
tion and the bulk viscosity arc proportional to sin Atp. 
(Note that the departure from the thermal equilibrium 
is negligible because it is restored by strong forces on 
much shorter time scales.) 

It should be clear that the instantaneous flavor compo- 
sition in oscillating quark matter and the rate difference 
of the forward and backward weak processes in Fig. Q] arc 
related to each other. The difference of the rates changes 
the composition, while the composition in turn influences 
the difference of rates. The corresponding dynamics can 
be conveniently described in terms of the time dependent 
deviations of the chemical potentials from their equilib- 
rium values. 

In (5 equilibrium, the chemical potentials of the three 
lightest quarks are related as follows: fi s = fid = fi u + Me- 
Here fi u , fid and fi s are the chemical potentials of up, 
down and strange quarks, while fi e is the electron chem- 
ical potential. In pulsating matter, the instantaneous 
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departure from equilibrium is described by the following 
two independent parameters: 



Sfii = fx s - Hd = 8fi s - Sfi d , 
S/J.2 = Ms - M« - Me = SfJ, a - 5fJ, v 



(5a) 
(5b) 



where 5/ii denotes the deviation of the chemical potential 
fii from its equilibrium value. (Note that 6/13 = \id~ \i u ~ 
fx e =8^2 — 8^1 is not independent.) When Sfii are non- 
zero, the corresponding pairs of forward and backward 
weak processes in Fig. Q] have different rates. To leading 
order, the rate differences are linear in 5/ii, 



r 



r 



(a) 
(e) ' 



- (6) 
"(d) 

'(/) 



-A 3 (5{i 2 - S^i) . 



(6a) 
(6b) 
(6c) 



The corresponding A-rates have been calculated for the 
normal phase [13, HH, [U [42| as well as several color su- 
perconducting phases of quark matter [H, H3, 5(J ■ The 
results for the normal phase, in particular, read 
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(7a) 
(7b) 
(7c) 



These will be used below as a benchmark for the rates in 
spin-one color-superconducting phases. 

The semi-leptonic rate A3 is determined by the Urea 
processes u + e~ — > d + v e and d — > u + e~ + v ei shown 
in diagrams (e) and (/) in Fig. [TJ It was calculated in 
Ref. [37| for four different spin-one color-superconducting 
phases of quark matter. The result has a form of the 



product of the rate in the normal phase A 3 
specific suppression factor, 



(0) 



A3 = A 



(0) 



2 -H, 
3 \T 



and a phase- 



(8) 



where is the spin-one color-superconducting gap pa- 
rameter, and H(cj)/T) is a suppression factor for the pro- 
cesses involving gapped quasiparticlcs. (The first term 
in square brackets is the contribution of ungapped quasi- 
particles.) When </>—>• 0, the suppression factor H(<j)/T) 
approaches 1 and the normal phase result is restored. A 
simple fit to the numerical data of Ref. [13 for H((p/T) 
is presented in Appendix [A] 

Because of similar kinematics and phase space con- 
straints for the other pair of semi-leptonic processes, 
u + e~ — > s + u e and s u + e~ + v e , shown in diagrams 
(c) and (d) in Fig. [I] the dependence of the rate A2 on the 
color-superconducting gap should take the same form as 
A 3 in Eq. (JS), i.e., 



A 2 — A 2 



(0) 



3 3 \T 



(9) 









FIG. 1: Diagrammatic representation of the /3 processes that 
contribute to the bulk viscosity of dense quark matter. 



In contrast, the rate Ai is determined by the nonleptonic 
processes u + d — > s + u and s + u — > u + d, see diagrams 
(a) and (b) in Fig. [TJ which have a qualitatively different 
kinematics. In spin-one color-superconducting phases of 
quark matter, this was recently calculated in Ref. [40j |. 
The numerical result can be conveniently summarized by 
the following expression: 



A, = Af) 



Af+(l-Af)H{ r. 



(10) 



where, in addition to the suppression factor H(<j>/T), we 
also introduced a constant Af, which determines a rela- 
tive contribution of the ungapped quasiparticles to the 
corresponding rate. In the four spin-one phases stud- 
ied in Ref. [4fJ, the constant takes the following val- 
ues: N A = Af polar = 1/9, A/" planar « 0.0393, and 
A/" CSL = 928/27027 0.0343. A simple fit to the nu- 
merical data for H(<f>/T) is given in Appendix [Bl 

When the rates ©, © and (fTU)) arc known, the calcu- 
lation of the instantaneous pressure and, thus, the bulk 
viscosity from Eq. is straightforward [Hj]. Here we 
quote only the final expression for the viscosity, 



C = Ci + C 2 + C 3 , 



111) 
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where 

Ci 



n a 2 a 3 



aia 2 a 3 Ci 



+ (ai + Q' 2 + a 3 ) (Aid - A 2 Ci) 
n ai«3 



(12a) 



w 5? + 52 



0:10203(72 

(01 + o 2 + 03) [(-42 - B 2 ) Ci - ^ 2 C 2 ] Z (12b) 
aio 2 a 3 (Ci - C 2 ) 2 



Cs = - 



n ai« 2 



w Si + S2 



+ (ai + Q'2 + a 3 ) (BiC 2 - B 2 dY 



(12c) 



and 

Si = 
32 = 



—010-203 + (Q'i + Q'2 + Q'3) (BiA 2 - AiB 2 ) , (13a) 
aio 2 (Bi - B 2 ) + aia 3 (A 2 - B 2 ) + a 2 a3j4i.(13b) 



Here a, = nu/\i (i = 1,2) and n is the baryon den- 
sity of quark matter. The quantities Aj. Bi and Ci are 
susceptibility-like functions, see Ref. [34[ for the defini- 
tion. To leading order in (fr/ni, they are the same as in 
the normal phase. 

For comparison, let us also note that the bulk viscosity 
in the limit of the vanishing scmi-lcptonic rates reads 



u) aj+ A{ 



(14) 



III. 



NUMERICAL RESULTS FOR BULK 
VISCOSITY 



In our calculation of the bulk viscosity in spin-one 
color-superconducting quark matter below, we choose the 
same two representative sets of model parameters as in 
Ref. M- 



Set A 


Set B 


n = 5/o 


n = 10po 


m s = 300 MeV 


m s = 140 MeV 


a 3 = 0.2 


a s = 0.1 



In both cases, the masses of light quarks are the same: 
m u = 5 MeV and rrid = 9 MeV. In accordance with gen- 
eral expectations, the values of the strange quark mass 
m s and the strong coupling constant a s should be larger 
(smaller) in the case of lower (higher) density. This qual- 
itative property is reflected in the model parameters in 
Set A (Set B). The values of all chemical potentials as 
well as the coefficient functions Ai, Bi and Ci for each 
set of parameters are quoted in Table HI 

It may be appropriate to briefly comment about the 
choice of the strong coupling constant a s in the model at 
hand. The values of a s in both sets of parameters may 
seem abnormally small. Indeed, the running coupling in 



QCD is about 0.12 at the scale of Mz (mass of Z boson) 
and about 0.32 at \/3 GeV [43| . However, here we use the 
model parameter a s only in order to capture several qual- 
itative (Fermi liquid) effects in quark matter. Its nonzero 
value allows (i) to avoid the underestimation of the rate 
of semi-leptonic processes due to a limited phase space 
[24I ] and (ii) to mimic the modification of the quark equa- 
tion of state due to strong interactions, see Ref. [34| for 
details. The naive extension of the corresponding leading 
order corrections to the regime of strong coupling is prob- 
lematic. Not only this would imply the use of the pertur- 
bative results beyond the range of their validity, but this 
would also lead to very large and seemingly unphysical 
effects on the equation of state, used to determine the 
susceptibility functions Ai, Bi and Cj. (Notably, if the 
equation of state is kept unchanged, the increase of a s in 
the A3 -rate, even by an order of magnitude, has little 
effect on the viscosity.) This dilemma could be resolved 
by properly accounting the non-perturbative dynamics 
of QCD. At present, however, such a task seems insur- 
mountable at the low energy scales relevant for neutron 
stars. For the purposes of this study, therefore, we treat 
a s as a small independent parameter that captures only 
some qualitative properties of quark matter. 

Here the critical temperature of the spin-one color- 
superconducting phase transition is assumed to be T c = 
2 MeV. This may be a somewhat high, but still reason- 
able value for T c . Indeed, in QCD the spin-one gap is 
estimated to be about two orders of magnitude smaller 
than the spin-zero gap [20[-22l. and the latter is natu- 
rally of order 100 MeV [l4Ml9l |. Even higher values of the 
spin-one gap have been reported in Ref. 



441. The effect 



of varying the critical temperature is easy to understand 
and will be briefly discussed below. As in Ref. [33], we 
use the following model temperature dependence of the 
gap parameter: 



<j>(T)=<j> dl- 



for T < T c 



(15) 



with 0o being the value of the gap parameter at T = 0. 
Note that the ratio T c /cf>o depends on the choice of the 
phase [l^ . The approximate values of this ratio are 0.8 
(CSL), 0.66 (planar), 0.49 (polar), and 0.81 (4-phase). 

For model parameters in Set A, the numerical results 
are presented in Fig. [5] As we can see, the value of T c 
determines the point where the bulk viscosity starts to 
deviate from the benchmark result in the normal phase 
(shown by the gray solid line). The upper panels show 
the dependence of the bulk viscosity £ on temperature 
for two representative values of the oscillation frequency, 
t -1 = 10 Hz and r _1 = 1000 Hz. The lower panels in 
the same figure show the temperature dependence of the 
ratio C/Cnon, where £ is the bulk viscosity that takes into 
account all weak processes, while C non is an approximate 
result, see Eq. (I14[) . in which only the nonleptonic pro- 
cesses are included and the semi-leptonic processes are 
not. When the ratio C/Cnon is substantially larger than 
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TABLE I: Two sets of parameters used in the calculation of the bulk viscosity. 



model 


Me [MeV] 


fin [MeV] 


[i d = ^ [MeV] 


Ai [MeV] 


A 2 [MeV] 


B 1 [MeV] 


B 2 [MeV] 


Ci [MeV] 


C* 2 [MeV] 


Set A 


39.139 


402.463 


441.602 


239.432 


127.937 


111.386 


-3.726 x 10 4 


-60.463 


-60.460 


Set B 


7.396 


495.275 


502.671 


324.556 


164.288 


160.268 


-2.080 x 10 6 


-10.692 


-10.709 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 



T [MeV] 




0.01 0.05 0.10 0.50 1.00 5.00 10.00 

T [MeV] 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 



T [MeV] 



FIG. 2: (color online). Temperature dependence of bulk viscosity £ and the ratio C/Cnon for model parameters in Set A and the 
spin-one color-superconducting critical temperature T c = 2 MeV. The results for two fixed frequencies of the density oscillations 
are shown. 



1, it is an indication that the semi-lcptonic processes play 
an important role and, thus, cannot be neglected. 



Compared to the normal phase result, the main fea- 
tures of the temperature dependences in spin-one color- 
superconducting phases (see the upper panels in Fig. [2]) 
are (i) a smoothed shape of the semi-leptonic "hump" 
and (ii) an overall enhancement of the bulk viscosity due 
to color superconductivity for a substantial range of tem- 
peratures below T c . 



As in the case of the normal phase, the semi-leptonic 
processes are responsible for an increase ( "hump" ) of the 
bulk viscosity in a region of temperatures around Thump, 



where 



^ P A) - 2.1 MeV ' 



T^± B) ~ 1.4 MeV 



" hump 



1 ms 



1/4 



1/4 



(16a) 
(16b) 



are the approximate positions of the peak of the hump 
in the normal phase in the case of the model parameters 
in Set A and Set B, respectively. In order to derive these 
results, we used an approximate expression for the bulk 
viscosity in Eq. (24) of Ref. (34[, which is valid when the 
nonleptonic rate is infinitely large while the semi-leptonic 
rates are finite. The maximum of that expression corre- 
sponds to A 2 + A 3 = nujAi/(BiA 2 — B 2 Ai), whose so- 
lution determines an approximate value for Thump- Two 
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FIG. 3: (color online). The contour plot of the bulk viscosity enhancement factor due to spin-one color superconductivity. The 
results are for the CSL phase in a model with the parameters in Set A. The results for two frequencies of the density oscillations 
are shown: r _1 = 10 Hz (left panel) and r _1 = 1000 Hz (right panel). The ratio CcSL/Cnormal is larger than 1 in the colored 
regions and is equal to or less than 1 in the white region. The contours are labeled by the corresponding values of the ratio 

CcSL/Cnormal - 



remarks are in order here: (i) the scaling law Thump <x 
1/t 1 / 4 follows from the power-law temperature depen- 
dence of the semi-leptonic rates A 2 , A3 oc T 4 and (ii) the 
overall value in Eq. ((16)) is slightly corrected to match the 
actual numerical results in the case of a finite nonleptonic 
rate. 

When T c > Ti lump the semi-leptonic hump is partially 
washed out by the presence of color superconductivity. 
This is most clearly seen from the ratio of the bulk vis- 
cosities C/Cnon in the lower panels in Fig. [2] While the 
inclusion of the semi-leptonic processes leads to an in- 
crease of the viscosity, the effect is not as large as in the 
normal phase. Of course, this conclusion is sensitive to 
the choice of the color-superconducting critical temper- 
ature T c . In general, two qualitatively different regimes 
can be realized. When T c < Thump; the hump occurs in 
the normal phase and, therefore, its shape is almost un- 
affected by color superconductivity. In the opposite case, 
T c > Ihump, the effect is present and gets stronger as T c 
increases relative to Thump- 

Now, let us turn to an overall enhancement of the bulk 
viscosity due to color superconductivity below T c . This 
is observed almost for the whole range of temperatures 
To, max ^ T < T c , where To, max is the temperature at 
which the bulk viscosity of the normal phase has a global 
maximum. The value of T^max can be easily estimated by 
considering an approximate expression for the bulk vis- 
cosity (fTf|) when only the nonleptonic processes are taken 
into account. The maximum of Eq. (|14p corresponds to 
a.\ = A±. After solving this for the temperature, we ob- 
tain 

«ax A) - 47 keV^, (17a) 
ifiL? - 41 keV^, (17b) 



where r is the period of oscillations measured in millisec- 
onds. Notably, the location of the maximum is almost the 
same for both sets of model parameters. Because of the 
superconductivity, the location of the maximum is shifted 
to a higher temperature, T^max — To,max/vW, where Af 
is the same parameter that appears in Eq. (|10[) . Taking 
the shift of the maximum into account, we find that the 
enhancement relative to the normal phase is observed for 
T±<T <T C with T* ~ To.max/A/" 1 / 4 being the point be- 
tween TJ^max and T0 jmax , at which the bulk viscosities for 
the normal and superconducting phases cross. At lower 
temperatures, T < T±, the effect of color superconduc- 
tivity is opposite: it reduces the bulk viscosity. 

The range of temperatures, in which the bulk viscos- 
ity increases relative to the normal phase of quark mat- 
ter, depends on the value of the critical temperature T c 
and the frequency of oscillations. While the actual en- 
hancement of the viscosity also depends on the specific 
pattern of spin-one pairing, the qualitative features in 
all four phases studied here are similar. As an example, 
let us consider the CSL phase in more detail. In Fig. [3J 
we show the contour plot for the bulk viscosity enhance- 
ment factor due to color superconductivity. The ratio 
CcSL/Cnormal is larger than 1 only in the colored regions 
in Fig. |31 In white regions, it is either 1 (when T > T c ) 
or less than 1 (otherwise). 

As evident from Fig. [JJ the enhancement of the bulk 
viscosity by spin-one color superconductivity occurs in a 
rather wide range of temperatures, especially when the 
frequency of density oscillations is not too large and the 
value of T c is not too small. At r _1 = 10 Hz, for ex- 
ample, it extends over an order of magnitude or more in 
temperature, provided T c > 100 keV. At t" 1 = 1000 Hz, 
in contrast, an order of magnitude or wider temperature 
range for the enhancement is seen only if T c > 1 McV. 
(It should be noted that, in the case r _1 = 10 Hz shown 
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FIG. 4: (color online). Temperature dependence of bulk viscosity £ and the ratio C/Cnon for model parameters in Set B and the 
spin-one color superconducting critical temperature T c — 2 MeV. The results for two fixed frequencies of the density oscillations 
are shown. 



in the left panel in Fig. [3J the ratio CcSL/Cnormai is truly 
less than 1 in a small white region just below the T = T c 
line. This "abnormality" is due to a subtle interplay be- 
tween the semi-leptonic and nonleptonic processes when 
the value of T c is fine-tuned to be near Thump-) 

By ignoring the subtle complications due to the semi- 
leptonic hump around Thump, we find that the enhance- 
ment of the bulk viscosity in the window of temperatures 
T* < T < T c (as well as the suppression at lower tem- 
peratures, T < T*) is primarily due to the reduction of 
the nonleptonic rate Ai in color-superconducting phases. 
At temperatures below T c , when all gapped quasipar- 
ticles effectively cease to contribute, the corresponding 
reduction factor for the rate is approximately given by 
the value of Af. This means that the enhancement factor 
for the viscosity approaches its inverse value, A/" -1 . By 
making use of the numerical results for A/", we find that 
the enhancement factor for the bulk viscosity reaches up 
to about 9 in the A- and polar phases, 25 in the planar 
phase and 29 in the CSL phase. (The suppression factors 
at T < Tt approach the same values.) In the region of 
the hump, of course, the behavior is more complicated, 
but the overall effect of superconductivity is still mainly 
to increase the bulk viscosity. 



The numerical results in the case of the model parame- 
ters in Set B are shown in Fig. [4] The qualitative features 
are similar to those obtained for Set A. However, the ef- 
fect of the semi-leptonic processes is less pronounced: the 
corresponding hump is almost non-existent and the ra- 
tio C/Cnon does not much deviate from 1. At the same 
time, the effect of color superconductivity is very well 
pronounced. Compared to the normal phase result, an 
enhancement of the bulk viscosity by a factor of about 
A/"" 1 is seen in a relatively wide window of temperatures 
from T- to T c . 



IV. DISCUSSION 

In this study, we calculated the bulk viscosity in spin- 
one color-superconducting strange quark matter by care- 
fully taking into account the interplay between the non- 
leptonic and semi-leptonic week processes. 

As expected, the nonleptonic processes give the domi- 
nate contribution to the viscosity in a wide range of pa- 
rameters. Yet, as in the normal phase [34j], the semi- 
leptonic processes may also lead to a substantial cor- 
rection in a window of temperatures around Thump, see 



8 



Eq. (fT6|) . The value of Thump scales as 1/r 1 / 4 and hap- 
pens to be of order 1 MeV for millisecond pulsars. The 
size and the relative importance of the hump can be con- 
veniently measured by the ratio C/Cnon when it is notice- 
ably larger than 1. For millisecond pulsars, however, this 
ratio remains close to 1. The effect is more pronounced 
when the period is a few orders of magnitude longer. We 
also find that the corresponding hump in the temperature 
dependence of the bulk viscosity of color superconductors 
is partially washed out compared to the normal phase. 
The higher is T c relative to Turnip, the larger wash out of 
the hump is seen. At sufficiently low T c , i.e., T c < Thump, 
the hump occurs in the normal phase and, therefore, its 
shape is unaffected by color superconductivity. 

If the critical temperature of the spin-one color- 
superconducting phase transition T c is larger than To jmax , 
see Eq. p7|) . the main effect of color superconductivity is 
an overall increase of the bulk viscosity in a range of sub- 
critical temperatures, T+ < T < T c , sec Fig. [31 The cor- 
responding range of temperatures widens with increasing 
the value of T c and with decreasing the frequency of oscil- 
lations. The increase of the viscosity is primarily due to 
the suppression of the nonleptonic rate by color supercon- 
ductivity. At almost all temperatures below T c , the rate 
is dominated by the ungapped quasiparticles, whose rel- 
ative contribution is scaled by the factor Af with respect 
to the normal phase (note that J\f < 1). It is the inverse 
value A/" -1 that determines the maximal enhancement of 
the bulk viscosity at subcritical temperatures. The corre- 
sponding enhancement factor is equal to 9 in the A- and 
polar phases, about 25 in the planar phase and about 
29 in the CSL phase. (At temperatures below T+, color 
superconductivity leads to a suppression of the bulk vis- 
cosity, and the maximal suppression will approach the 
same value of A/"" 1 .) 

In relation to this result, it might be appropriate to 
note that a similar enhancement mechanism was previ- 
ously observed for spin-zero color superconductors [36j . 
A special feature of spin-one color superconductivity is 
that the maximum enhancement factor can be much 
larger. 

In our analysis, we utilized the same spin-one pairing 
pattern as in Refs. [20U22I ] . In the case of zero quark 
masses, the main signature of the corresponding phases 
is the presence of ungapped quasiparticles. When quarks 
have small masses, the gaps of the corresponding modes 
are of order tpm/ /i. These may be still too small to signif- 
icantly affect our main results. However, if the spin-one 
gaps are larger, as some studies suggest j44j, the sup- 
pression of the nonleptonic rates and, therefore, the en- 
hancement of the bulk viscosity in color superconducting 
matter may turn out to be even stronger. 

In application to compact stars, we may speculate that 
the transition to a spin-one color superconducting phase 
in a stellar core can have a stabilizing effect against the 
r- modes driven by the gravitational radiation [25J . If the 
critical temperature of the corresponding phase transi- 
tion is on the order of or above 1 MeV, the correspond- 



ing dynamics can affect even relatively young stars. The 
study of the actual quantitative effect that spin-one color 
superconductivity has on the reduction of the instability 
window in the pulsar frequency and temperature plane 
can be done along the lines of Refs. |45T - [47j |. However, 
this problem is beyond the scope of the present paper. 
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Appendix A: A-rates of semi-leptonic (Urea) 
processes 

The rates of the semi-leptonic processes in spin-one 
color superconducting quark matter were calculated in 
Ref. [I?}. The general expression for the rate takes the 
following form: 



A; = A 



(0) 



5 + 



for 



2,3, 



(Al) 



where A,^ is the corresponding rate in the normal phase 
of quark matter and H {4>/T) is a phase-specific suppres- 
sion factor. By construction, it satisfies the constraint 
-ff(O) = 1, which corresponds to the case of the normal 
phase. We used the numerical data of Ref. [33] to obtain 
the following fits for the suppression factors as functions 
of the dimensionless ratio tp = 4>/T in the four spin-one 
color superconducting phases of quark matter: 



H A (tp) 



a\tp A + blip 3 + ciLp 2 + d\ 



(A2) 



Lp 5 + dip 3 + fnp 2 + di 

where a x = 1.069, h = -0.2187, a = 3.666, d x = 21.50, 
e x = 1.333 and /1 = 9.349, 

a 2 tp 3 + b 2 tp 2 + c 2 



if polar (tp) 



d 2 tp 4 + e 2 tp 3 + f 2 ip 2 + c 2 ' 



(A3) 



where a 2 — n, b 2 = 21.94, c 2 
11.20 and h = 214.0, 



^planar ^ 



a 3 tp 



1386, d 2 = 6.994, e 2 

5 + b 3 tp 3 + c 3 tp 2 + d 3 (l + tp) 



<p 



tp 3 + e 3 tp 2 + d 3 



(A4) 

where a 3 = 0.917, b 3 = 0.456, c 3 = 11.69, d 3 = 34.0 and 
e 3 =4.221, 



H 



CSL 



q 4 ^ 4 + bjtp 3 + dtp 2 + dj(l + V2tp) _^ 2 - lf> 



tp 3 + eitp 2 + d 4 



(A5) 
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where a 4 = 1.034, 64 = 1.001, c 4 = 9.735, ft = 13.81 and 
e 4 = 1.684. 



Appendix B: A-rates of nonleptonic processes 

The A-rate of the nonleptonic processes in spin-one 
color superconducting quark matter was calculated in 
Ref. [i^. The general expression for the rate takes the 
following form: 



Ai = \\ 



(0) 



Af + (l-Af)H{r 



(Bl) 



where xf^ is the corresponding rate in the normal phase 
of quark matter, J\f is a constant that determines the 
relative contribution of ungapped quasiparticles to the 
rate, and H (4>/T) is a phase-specific suppression factor 
due to gapped quasiparticles. The normal phase corre- 
sponds to (j) = 0, in which case there is no suppression 
and H(0) = 1. The value of N for each phase reads 



N A 

^"planar 

A/" CSL 



1 

9' 
1 

9' 

0.0393, 
928 
27027' 



(B2a) 

(B2b) 
(B2c) 
(B2d) 



For this study we used the numerical data of Ref. [40J 
to obtain the following fits for the suppression factors as 
functions of the dimcnsionless ratio ip = <fi/T: 



QlV3 2 + ft 

ip 3 + lip 2 + Pi 



(B3) 



where a x = 0.1247, ft = 12.60 and 71 = 5.042, 



a 2 tp 2 + ft 

(ft + 72<ft + ft ' 



(B4) 



where a 2 = 0.0271, ft = 65.45 and 72 = 13.35, 



jyplanar ^ = 



C*3<P 4 + fty 3 + 73V? 2 + ft(l + ft 

<ft + ft 



(B5) 

where a 3 = 0.0717, ft = -0.2663, 73 = 1.108 and ft = 
4.561, 



H CSL (ft) 



c^tft + ftjft + 74 (1 + y/2<p) 
ft + ft ft + 74 



-V2tp 



(B6) 



where a 4 = 0.6981, ft = -2.045, 74 = 4.482 and ft 
-1.217. 
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